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ABSTRACT 

Balmer-dominated shocks in supernova remnants (SNRs) produce strong hy- 
drogen lines with a two-component profile composed of a narrow contribution 
from cold upstream hydrogen atoms, and a broad contribution from hydrogen 
atoms that have undergone charge transfer reactions with hot protons. Obser- 
vations of emission lines from edge-wise shocks in SNRs can constrain the gas 
velocity and collisionless electron heating at the shock front. Downstream hydro- 
gen atoms engage in charge transfer, excitation and ionization reactions, defining 
an interaction region called the shock transition zone. The properties of hot hy- 
drogen atoms produced by charge transfers (called broad neutrals) are critical for 
accurately calculating the structure and radiation from the shock transition zone. 
This paper is the third in a series describing the kinetic, fluid and emission prop- 
erties of Balmer-dominated shocks, and is the first to properly treat the effect 
of broad neutral kinetics on shock transition zone structure. We use our mod- 
els to extract shock parameters from observations of Balmer-dominated SNRs. 
We find that inferred shock velocities and electron temperatures are lower than 
those of previous calculations by < 10% for v s < 1500 km s -1 , and by 10 — 30% 
for v s > 1500 km s -1 . This effect is primarily due to the fact that excitation by 
proton collisions and charge transfer to excited levels favor the high speed part of 
the neutral hydrogen velocity distribution. Our results have a strong dependence 
on the ratio of electron to proton temperatures, (3 = T e /T p , which allows us to 
construct a relation (3{y s ) between the temperat ure ratio and shock veloci ty. We 



compare our calculations to previous results by iGhavamian et al.l (l2007bl ) 
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Introduction 



Balmer-dominated shocks in supernova remnants (SNRs) encounter upstream gas con- 
taining a substantial fraction of neutral hydrogen atoms. Emission from the shock is char- 
acterized by strong Balmer and Lyman lines with two-component profiles, which consist of 
a narrow contribution from direct excitation of the hydrogen atoms entering the shock, and 
a broad contribution from excitation of "broad neutrals," hydrogen ato ms that have been pro- 
duced by charge transfer rea ctions with protons in the downstream gas (j Chevalier fc Raymond 
1978 ; Chevalier et al] 1980). Such s hocks are seen in many SNRs, including par ts of the 



Cygnus Loop (IGhavamian et al. 



19871 . KCW87; iSmith et al.lll991l 



200ll . here after G01), Tycho and Kepler's remnant s (IKirshner et al 



, S91; G01 



2007bl . G07b), SN 1006 (KCW87; S91 ; IGhavamian et al 



Fesen et al. 



989, F89), RCW 86 (G01: IGhavamian et al 



2002 



, G02) , an d sev eral remnants lo- 



cated in the LMC flTuohv et alll982 . T82; S91; IGhavamian et al.ll2003l . G03; IGhavamian et al 
2007al . G07a). 



Profiles of broad emission lines from edge-wise observations of shocks in SNRs can be 
used to infer shock velocities. These calculations can b e compared to analyses of shock front 



prop er motion to compute distances to the objects (jChevalier et al.l Il980l ; IKirshner et al. 
19871 ). Observations which resolve the spatial profile of the combined Ha emissio n can 



constrain the neutral fraction and density of the upstream gas ([Raymond et al.l 120071 ) . 



In addition to diagnosing parameters of SNRs, Balmer-dominated shocks provide an 
important probe into the plasma physics of collisionless, non-relativistic shocks. Models of 
the broad component width and integrated broad-to-narrow intensity ratio can be used to 
derive the ratio of electron to proton temperature as a function of shock velocity. Such a 
relation can provide insight into the physical mechanisms at work in the collisionless plasma. 

This is the third in a series of papers in vestigating the hydrodyn amics, kinetics, and line 
emission from Balmer-dominated shocks. In lHeng fc McCrayl (120071 . Paper 1), we calculated 
velocity distribution functions for the broad neutrals and computed t h e rati o of broad-to- 
narrow line emission as a function of shock velocity. In iHeng et al.l (120071 . Paper 2) we 
calculated the density structure of the shock transition zone, where hydrogen atoms passing 
through the shock front undergo charge transfer reactions, emit radiation, and become fully 
ionized. 

The approximations employed in Paper 2 limit the validity of the results to shocks 
entering the upstream gas with velocities < 3000 km s -1 , due to our treatment of the broad 
neutrals as a single fluid with the same bulk velocity as the ions (the "restricted three- 
component model" of Paper 2, §3). In fact, the broad neutrals in the shock transition 
zone are not a fluid, but have distinct anisotropic distribution functions depending on the 
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number of charge transfer reactions they engage in (Paper 1). In the present paper, we 
treat charged species in the shocked gas as fluids, and describe the hydrogen atoms with 
appropriate kinetic distribution functions. In addition, we provide an improved calculation 
of the broad line velocity profile for v s > 2000 km s _1 . Using this methodology, we calculate 
the structure of the shock transition zone more accurately than in Paper 2, and characterize 
hydrogen line emission from shocks with 300 < v s < 10, 000 km s _1 . 

We find that we can self-consistently determine shock parameters for most Balmer- 
dominated SNRs. Our results yield lower values of the inferred shock velocity and proton- 
electron temperature equilibration than those derived using previous models. We compute 
the dependence of the elec tron temperature on shock velocity and compare it to the results of 
Ghavamian et al. J2007bh . We note that our calculations are unable to fit the observations 



for several SNRs; in these cases, our basic model must be augmented by new physics to 
account for the data. 

In $2] of this paper, we describe our physical model. In $3j we display the equations 
employed to calculate the structure of the shock transition zone and describe our numerical 
solution. In §HJ we compute the spatial emission profile and hydrogen line spectra from 
the shock transition zone. In §0 we analyze our results and describe how observations of 
Balmer-dominated SNRs should be interpreted in light of our new calculations. Finally, in 
§6j we discuss implications of our results for collisionless electron heating, explore limitations 
of our model, and identify areas for future research. 



2. Physical Model 

We consider a shock with velocity 300 < v s < 10, 000 km s -1 traveling through the 
ISM, which consists of cold, partially neutral hydrogen and neutral helium. The pre-shock 
fraction of helium, relative to hydrogen, is denoted fn e - If the total upstream density of 
protons and hydrogen atoms is n , then the fraction of pre-shock protons is defined to be 
f P = n p /n . 

In the frame of the shock, the cold upstream neutrals and charged particles flow down- 
stream with uniform velocity v s . At the shock, we assume that the protons are heated in a 
thermal distribution to T p ~ 3 j^ s ~ 10 7 K for v s ~ 10 s cm s _1 . At the high post-shock 
temperatures in Balmer-dominated remnants, emission from hydrogen exc itation is much 



stron ger than forbidden transitions in metals seen in radiative shocks (e.g., IShull &: McKee 



19791 ). The length scale for thermalization of protons is set by the proton cyclotron gyrora- 

gyro 



dius, /g yr o ~ 10 8 cm at B ~ 10 4 G, which is much less than the length scale for the shock 
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transition zone, / zone > l/(^o°") = 10 14 cm (set by the mean-free paths for excitation, charge 
transfer and ionization reactions, and the total upstream density). 

Without any energy transfer between shocked protons and electrons, the electron tem- 
perature is expected to be a factor of m e /m p less than that of the protons. Since the 
collisional timescale for temperature equilibration between elect rons and protons is in some 



cases longer than the age of the remnant (e.g., iSpitzerl Il962l ). we assume that the elec- 
trons remain at a fixed temperature with respect to the protons in the downstream material 
(note, however, that intra-species and inter-ion collisional timescales can be much shorter) [j] 
However, several authors have argued that electrostatic instabilities at the shock front can 
increase the electron tempera t ure, with estimates rang ing from T e ~ 0.1T p to T e = T p (see 
Cargill fc Papadopoulod Il988t iGhavamian et al.ll2007bl . and the references therein). Since 
the physics of non-relativistic, collisionless shocks is still poorly understood, we parametrize 
the ratio of electron to proton temperatures using the definition 

T 



T 



1 



One of the goals of our work is to use observations of Balmer-dominated remnants to con- 
strain the value of (3 in the shock transition zone as a function of shock velocity. 

In the uniform velocity upstream gas, we assume that negligible interactions take place 
between neutral and charged species. The cold neutral hydrogen atoms passing through the 
shock are not affected by the discontinuity. In contrast, the bulk velocities of the downstream 
thermal protons and electrons become approximately v s /4, at which point ionization and 
charge transfer reactions between cold atoms and hot protons occur (see §5.ip . However, if a 
significant amount of the energy dissipated in the shock is used to produce cosmic rays, the 
ion compression ratio increases, and their bulk velocity decreases. Furthermore, the presence 
of a cosmic ray precursor can broaden the distribution of upstream neutrals. We discuss the 
implications of these physical effects in §H1 

Downstream charge transfer reactions between cold neutrals and hot protons produce 
a new population of hot atoms, referred to as broad neutrals. Initially, such an interaction 
will produce a cold proton, which is re-energized through gyro motion around the magnetic 
field and intra-species Coloumb collisions with the hot proton population. We assume that 
su ch protons rapid l y equ ilibrate back into the thermal pool. However, recent calculations 
by iRaymond et al.l (120081 ) indicate that interactions within the shock front may produce a 
distinct "pickup" ion population analogous to that in the solar wind. 



1 Somc collisional heating of electrons by protons will occur. We estimate that such an effect can lower 
the derived value of (3 by f — 2% 
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Ha emission is produced by excitation reactions and charge transfers to excited states 
in the shock transition zone. The narrow and broad components of the line are emitted by 
the cold hydrogen and broad neutrals, respectively The spatial and velocity distributions 
of the broad neutrals are critical for calculating the correct structure and emission from the 
transition zone. In Paper 2, we assumed that the broad neutrals could be characterized 
as a single fluid, with velocity and temperature equal to that of the ions. This treatment 
is flawed because, unlike the ionic species, the broad neutrals have negligible interactions 
among themselves. Therefore, each time a neutral engages in charge exchange, it becomes 
part of a new distribution with distinct velocity and temperature (see Paper 1) and is not 
subsumed into the original population. This fact has not, until now, been properly taken 
into account in the literature. 

The detailed kinetic properties of the particle distributions and their reaction rates 
were calculated in Paper 1. We use these methods, with updated values for the atomic cross 
sections (see Appendix |A]) , in our fluid calculation of the shock transition zone structure. 
In this picture, the broad neutrals form an infinite set of atomic populations with distinct 
reaction rates. In practice, as noted in Paper 1, the velocity distribution function of atoms 
experiencing many charge transfers rapidly converges to that of the protons, usually after 
only three or more interactions. Thus, only three distinct broad components interact with 
the other species. The bulk velocities and temperatures of the broad components are written: 



where v k and Tk are the velocity and temperature of the broad neutral population after 
k charge transfers. The coefficients F v k(v s , (3, f p , fue) and FTk{v s , (3, f p , fne) are calculated 
using the methods described in Paper 1 and, in principle, are also functions of v p and T p . 
However, we show in §5.1l that the proton velocity and temperature do not vary significantly 
enough in the shock transition zone to affect the values of these coefficients. For k > 3, F v k 
and Fxk are taken to be unity. 

Ionization of helium in the shock transition zone produces singly-ionized He + and alpha 
particles. We assume that the helium ions and protons are coupled through rapid inter- 
ion Coloumb collisions. In this case, all of the ions share the same bulk velocity (vi) and 
temperature (Tj). Figure [1] schematically displays the density variation of the different 
particle species in the shock transition zone. 



"Vk F v kV pi 

Tk = F Tk T p , 



(2) 
(3) 
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3. Spatial Structure of the Transition Zone 

3.1. Basic Equations 

We employ a plane-parallel coordinate system in the frame of the shock, in which the 
shock front is at z = 0. The density structure of the transition zone is determined by 
conservation of mass flux: 

d a 

— (n H v s ) = -n H ^ n s R iH ,s - n H n p R TQ ,p, (4) 

s=e 

d a 

— (niui) = n H n p R ToiP - n x y~] n s R iBtS - n x n p R Tt%n (5) 



dz 

d_ 

dz 



(n 2 v 2 ) = n 1 n p R TiP -n2^2n s R iB>s -n 2 n p R T>p , (6) 



d a 

(n k v k ) = n k _in p R T , p - n k ^ n s R iBtS - n k n p R TjP , (7) 



dz 



dz 

d 

dz 



d a 

(n H eV s ) = -nHe^2risRiHe,s, (8) 



(n He +Vi) = n He } j n s R iH e,s - n He + n s R iHe + ;S , (9) 



-e s=e 
a 



— (n a Vi) = n He +^2n s R iHe+tS , (10) 



dz 
Tz 



s=e 

N 



n H R iH . s + ^2n k R iB)S , (11) 



k=i 



where the rate coefficients R for the atomic interactions (in units of cm 3 s _1 ) are defined in 
Appendix |A] Charge conservation requires the electron density to obey the relation n e = 
n p + n He + + 2n a . The subscripts iH, iHe, and iHe + denote ionization of hydrogen, helium, 
and singly-ionized helium, while the index s runs over the charged particles participating 
in the reactions, s = {e,p, He + ,a} (electrons, protons, singly-ionized helium and alpha 
particles). The subscript T indicates charge transfer of cold hydrogen atoms with protons. 
As discussed in ^|2l the subscript k denotes the neutral hydrogen population which has 
engaged in k charge transfers, where the index runs from k = 1 to k = oo. While each broad 
population will have a separate rate coefficient for every atomic interaction, the values are 
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not sensitive to the details of the broad distribution functions. Therefore, we use one rate 
coefficient to describe the interactions of all the broad populations: the subscripts iB and T 
denote ionization and charge transfer, respectively. We list the references for our interaction 
cross sections and rate coefficients in Appendix |A] 

The system of differential equations is completed by expressions for conservation of 
momentum and energy flux: 



d 

dz 



m p n H v 2 s + ^ ( p k + m p n k vl) + Am p n He v 2 s + ^ (P s + m s n s v 2 ) 

k=l s=e 
1 °° / 1 \ 

-m p n H vl + ^2 ( p k + u k + 2 m p n k v l J V k+ 

k=l ^ ' 

2n He v 3 s + (p* 



0. 



(12) 



d_ 

dz 



U s + -m s n s v\ 



(13) 



where P = nksT is the pressure and U = ^nksT is the energy density. The ratios Pk/n^ and 
Uk/rik are set by the kinetic calculation of the broad neutral distribution functions through 
eq. dSJ), which assumes Tj = 3m p v 2 /lQ(l + /3)ks- Though the ion temperature varies by 
as much as 10% in the shock transition zone, this produces a negligible effect on the rate 
coefficients. Thus, the value of T& is proportional to the local ion temperature. Furthermore, 
our steady-state fluid calculation neglects kinetic evolution of the broad neutral distributions, 
which we expect to be a modest effect. 



We define a natural length scale for the transition zone, according to the formula: 

L = 



n R' 



(14) 



8 cm 3 s 1 . 



where R is a typical value for the interaction rate coefficents, taken to be 10 
We also define a set of dimensionless variables with \ = z/l z , r\ = n/n , u = v/v s , and 
e = ksT '/(m p v 2 ). As discussed above, broad neutrals which have engaged in three or more 
charge transfers can be described with the thermal proton distribution function. We therefore 
define r/ N = ^^3%- Summing over eqs. (J7J) with k > 3 and using these definitions, the 
formulas for the shock transition zone structure can be written: 



d 

dx 



d_ 



-1- (V2U2) 

dx 



s=e 

a 

s=e 
a 



-iB,s 



V2VpK 



T,p, 



(15) 



(16) 



(17) 
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d 

dx 



(VNUi 



d_ 

dx 



VHe 



d_ 

dx 



(VHe+Ui 



Tx {r]aUi 



d 

dx 



d 

dx 



rj2Tj p 7lT,p -VN/] VsTliB,s, 
s=e 

a 

iHe.s, 

s=e 

a a 

iHe,s — VHe+ iHe+,s 
s=e s=e 
a 

iHe+,si 

s=e 

a / N 



k=l 



d 

dx 



1\ a , 

( %efc + + ^ r]He + S ( Vses H — ~ 

k=l s=e V m P 

N a , 

+ { 5 VkUk + VkuT) + 4t] He + ^ ( ^VstsUi H S - 

k=l s=e \ m P 



where 1Z = R/R. 



(19) 
(20) 
(21) 
(22) 

(23) 
(24) 



3.2. Solution Method 

Eqs. (115p - (124l) are a set of ten coupled, non-linear ordinary differential equations, which 
are solved with the constraint of charge conservation. We compute eqs. (II 5 1) - (1221) directly 
in terms of the number density flux, denoted y = rju. The ion velocity and temperature, Ui 
and 6i, determine the broad velocity Uk and temperature through the relations (j^I)-©- 
Once y and u are known, the dimensionless number density rj = y/u can be substituted into 
the right hand side of the conservation equations. 

We solve eqs. (|23|) and (124")) subject to the boundary conditions y#(0) = 1 — f p , 
yne(0) = (1 - f p )f He , y P (0) = f P , and y k (0) = y H e+(0) = y a (0) = 0. The 
upstream neutral temperature is taken to be zero, while the protons and electrons are in 
thermal distributions with temperatures T u <C T p . Integrating these equations, substituting 
the initial conditions, and using the definition for y yields: 



N 



F ^y k + Yl —y° ^ + ^ + 4 ^ e - 1 - (i + p)f P e u - 4(1 - f p )f He \ Ui + 



v fc=l 



s=p 
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N p 

y^Vk + A/e + X>« ) = °' (25) 

v fe=l v ' k s=p 



N a \ / N 



Kkyk + J u i + ( 5 2 F ^ kVk + 5(3ye + 5 & ) e ' + 

fc=l s=p / \ fc=l s=P / 

+ 4t/ ffe - 1 - 5(1 + /?)/ p e u - 4(1 - Z,,)/^ = 0. (26) 

Note that the sums over charged species run from protons to alpha particles and do not 
include electrons, which appear separately in the equations (we neglect terms proportional 
to m e /m p ). The system (|25l)-(l26l) can be solved simultaneously to yield a quadratic equation 
in Ui with coefficients that are functions of x- The solution to the quadratic has two positive 
roots. Only one of the roots is less than unity; it represents the physical solution for the bulk 
ion velocity, which must be less than the shock velocity. We use a standard Runge-Kutta 
method to solve eqs. ( fT5l) - (l24l) . at each integration step using ( J25l) - (l26l) to compute the ion 
velocity which determines the particle densities. 

Figure [2] shows the density structure of the shock transition zone for v s = 1000 km s _1 , 
f p = 0.5, and [3 = 1. Using eq. f|T4l . the spatial coordinate has been converted from x to 
physical distance z behind the shock front. We use external density no = 1 cm -3 and helium 
fraction /# e = 0.1 for all calculations in this paper. The left and right panels of Figure [2] show 
the dimensionless densities for the neutral and charged species, respectively. The electron 
and proton densities have been scaled by a factor of 1/18. In the left panel, the solid curve 
shows a monotonic decrease in the density of the cold hydrogen atoms, which are removed by 
both charge transfer and ionization reactions. Charge transfer produces three populations 
of broad neutrals. At low velocities, charge transfer dominates over ionization, and many 
broad neutrals with k > 3 are produced. The dash-dotted curve shows the density of neutral 
helium, which is ionized farther downstream than hydrogen due to its smaller ionization 
rate coefficients. In the right panel, the solid curve shows singly-ionized helium, which is 
produced downstream the from neutral atoms, and is then ionized to yield a monotonically 
increasing population of alpha particles, shown by the dotted curve. The proton and electron 
densities are depicted by the short-dashed and long-dashed curves, which saturate when all 
of the neutral species have been depleted. The final electron density is slightly higher than 
that of the protons due to the presence of alpha particles. 



4. Line Emission 



Once the density structure of the transition zone is determined, we can compute the 
hydrogen line emission, including the spatial distribution and line profiles for the broad and 
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narrow components. We neglect collisional dexcitation and assume that every atom is excited 
from the ground state. Ha photons are produced by transitions from atomic levels 3s and 
3d to 2p, as well as from 3p to 2s. In the latter case, the atomic physics is complicated 
by a possible transition from 3p directly to Is, which results in a Lyman (3 photon. If the 
medium is optically thin to Ly/3 photons (Case A conditions), this possibility can be taken 
into account by proper weighting of the distinct angular momentum states in constructing 
excitation and change transfer cross sections for Ha emission: 



&Ha — C"3s + &3d + -£>3p,2sC3p, 



(27) 



where the factor B^ Pt 2 S ~ 0.12 is the fraction of transitions from 3p to 2s. If the medium is 
optically thick to Ly/3 emission (Case B conditions), re-absorption by ground state hydrogen 
effectively traps Ly/3 photons until they are re-emitted as Ha photons. In this case, all of 
the transitions eventually result in Ha emission, and we set B 3p ^ s ~ 1. 

Case A and B conditions represent the two extremes of media that are optically thin 
and thick to Ly/3 scattering. For the stationary atoms (i.e., cold hydrogen) which produce 
the narrow line, the optical depth to scattering of Ly/3 photons is, at line center, Tp ~ 
nnQ'Ln-np 1, with tlh ~ 0.5 cm" 3 , a 



10 14 cm 2 , and L nT1P ~ 



2 x 10 14 cm 



v see, e.g. 



Rybicki fc Lightmanlll979l ; ICoxil200ll ). The column of upstream neutrals will also line scatter 
Ly/3 emission, increasing the effective value of Tp. Thus, in conditions appropriate for many 
Balmer-dominated SNRs, partial scattering of Lyman ph otons will occur, produc ing results 
intermediate between Cases A and B for the narrow line (IGhavamian et al.ll200ll ). 



The conversion efficiency of narrow component Ly/3 to Ha was first computed by 
Chevalier et all (119801) . Subse q uent M onte Carlo calculations were performed by lLaming et al 



( 119961 ) and IGhavamian et al.l (120011 ). Briefly, models of the neutral hydrogen density and 



narrow component excitation rate are computed as a function of distance behind the shock. 
These are used to calculate the profile of excitations to the 3p level. Photons are emitted at 
frequencies distributed according to the pre-shock temperature in random directions. They 
are followed as they are absorbed and re-emitted as Ly/3 or Ha at a different locations un- 
til they escape. The conversion fraction depends on shock speed, electron-ion temperature 
ratio, and pre-shock ionization fraction, but not on total upstream density. We assume that 
the variation with /„ is modest, and that the temperature ratio decreases fr om (3 ~ 1 for 



< 500 km s 1 to (3 < 0.1 for v s > 2000 km s 1 , and use the calculations of lLaming et al. 



( 119961 ). in which 60 — 75% of 3p excitations result in Ly/3 photons. We fit the resulting 
conversion fraction as a function of shock velocity, yielding: 



B 3 p } 2s = 12vJ£ + (0.63 + 3.6 x 10 _ X, 



5), 



(28) 



where v s 5 is the shock velocity in units of 10 5 cm s 1 . 
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4.1. Spatial Emissivity Profile 

To determine the emission, we employ excitation rate coefficients calculated using the 
methods of Paper 1 (see also Appendix[A|). An emission line photon is produced when a cold 
hydrogen atom or broad neutral is excited by a charged particle, or undergoes charge transfer 
to an excited state. In the former case, the rates must be weighted by the probability of 
repeated excitation. We use the following equations to calculate the spatial emissivity profiles 
for the narrow (£ n ) and broad (£&) components: 

a 

= i y2 n sRAn,E ,s, (29) 



N ( 1 a ~ \ 

€b(z) = n H n p R ArhT * jP + ^ n k I n p R AnjT * iP + - — — ^2 n sR-An,E,s I 

k=l \ ^ s=e J 



(30) 



where the rate coefficients are labeled by the transition An = Ha, Lya, Lyf3, atomic inter- 
action, and particle type. The symbols E and Tq denote excitation and charge transfer 
to an excited state for cold hydrogen, while E and T* denote these interactions for broad 
neutrals. The probabilities Pe = Re / (Re + Ri + Rt ) and Pg = Re/ (Re + Ri + Rr) are 
calculated using the total reaction rates per atom or broad neutral for ionization, excitation, 
and charge transfer. These are calculated using the weighted sum of the rate coefficients, for 
example, i? 7 = Y,k=i n kY,^= e n sRiB,s- 

Figure [3] shows the spatial emissivity profiles for the narrow and broad components as a 
function of distance z behind the shock front. The thin vertical lines indicate the centroids 
z cn ^ of the emission components, calculated according to the formula: 

dz£, nib (z)- dz£ n!b (z)=0. (31) 

J z cn,b 

For the low value of f p shown in this figure, significant ionization must occur before there 
are enough charged particles to excite the narrow emission and engage in charge transfers 
to produce broad neutrals. Therefore, the intensity peaks for both components are shifted 
downstream from the shock front, and the centroid of the broad line emission is shifted 
farther downstream than that of the narrow line. 



4.2. Line Profiles 



The full width half maximum (FWHM) of the broad line can be related to the velocity 
and temperature equilibration of the shock. The line profile is a convolution of the broad 
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neutral and exciting species distribution functions with the cross sections for excitation and 
charge transfer to an excited state, projected along the line of sight to the observer. In a 
cylindrical coordinate system (r, 6, z) where the z axis is along the shock velocity direction, 
it is straightforward to project the line profile for observers oriented both edge-wise (along 
the r axis) and face-on (along the z axis) with respect to the shock front. Due to limb 
brightening, most observations will be selected for shocks viewed edge-wise or nearly so. 

We calculate hydrogen line profiles 4>foi an d 4>ew according to the formulas: 

4>Fo(Vz) 



4>Ew(Vr, Z) 



Fxs,z 
Fx.s,r 

where Av = |v — v'|, and the cross sections <Jx,Ha,s denote excitation and charge transfer re- 
actions for Ha emission to appropriately weighted angular momentum states, using eqs. (|27j) 
and (1281) . It should be noted that in eqs. (I32l - (l3~5l) . the quantities denoted by / are kinetic 
distribution functions, and should not be confused with the pre-shock ionization fraction 
f p . The edge-wise profile can be calculated as a function of position behind the shock. In 
practice, however, the observed lines are not spatially resolved. We therefore calculate the 
profiles as a function of v s and /3, spatially averaged over the shock transition zone. Figure H] 
shows examples of symmetric, edge-wise broad neutral velocity distributions, which are in- 
puts to eq. fl3"3"l) . Results are depicted in a reference frame where the average ion velocity 
is zero, for (3 = 0.01,0.1,0.5, 1, with v s = 1000 km s" 1 (left panel), and 7000 km s" 1 (right 
panel). For all cases, the pre-shock ionization fraction is set to f p = 0.5. 

For v s < 1500 km s _1 , excitation by electrons dominates over that by protons (including 
charge transfer into excited states). Even at moderate values of the electron temperature, 
(3 < 0.1, the thermal width of the electron distribution function is much larger than that 
of the broad neutral distribution, and the electrons typically have much higher velocities. 
Thus, for the majority of the range of integration, Av ~ v', and the electron contribution to 
eqs. f l32|) - fl33|) is approximately equal to the projected velocity distribution of the broad neu- 



n H npF T p, z (v s ) + 



dV 



^ n s f s (v')F E s,z{ v z> v ') + n pfp( v ')F Tp , z (v z , v') 



n H n p F Tpir (v s ) + 



n s f s (v')F E s,r{ V r> v ') + npfp(v')F T p,r{Vr, v') 

. s=e 

p2iv poo N 

/ d0 I dv r v r n k /fc(v) Av ax,H a ,s{^v), 
Jo Jo k=1 

p2it poo _^ _ 

/ d6 dv z n k f k (v) Av a x ,Ha,s {Av) , 
Jo J-oo k=1 



(32) 

(33) 
(34) 
(35) 



13 



trals, multiplied by the rate coefficient for excitatio n by electrons. This approximation has 



been used in previous studie s of the line profile (e.g.. IChevalier et al.lll980t iGhavamian et al. 



200ll ; lHeng fc McCrayll2007l ). Nevertheless, for shock velocities v s ~ 2000 km s , excitation 
rates by electrons and protons are comparable, with the proton contribution increasing for 
faster shocks. The cross sections for excitation by protons increase with the relative speed of 
the colliding particles. Since the broad neutrals and ions have comparable speeds, high speed 
neutrals are more likely to produce Ha photons. The effect on the line profile is somewhat 
mitigated by the integration over velocity space, but the observed line width is larger than 
the velocity width of the broad neutral distribution. This effect is relatively small at low 
velocity (< 10%) but is significant for high velocity shocks (see $6]). 

Figure [5] displays the FWHM of the broad, edge-wise Ha line profile as a function of 
v s , at several values of /3 — 0.01, 0.1, 0.5, 1. The pre-shock ionization is set at f p = 0.5. The 
FWHM increases monotonically with v s due to the increased temperature of the post-shock 
ion distributions (and hence broad neutral distributions) behind faster moving shocks. As 
f3 is increased, energy is transferred from the protons to heat the electrons, leading to lower 
proton temperatures and a smaller FWHM for the broad component. 

For fast shocks, in which proton excitation contributes substantially to eqs. fl32|) - fl33|) . we 
expect different transitions (e.g., Ha and Lya) to have different FWHM relations, since they 
employ distinct reaction cross sections in the calculation of <fi. This is a unique prediction 
made by our calculation, and may have consequences for studies of Lya emission from shocks 
with v s > 4000 km s -1 , in which the FWHM of the Lya line increases from 10-60% over 
that of the Ha line. 



5. Dependence on Shock Parameters 

Here we describe how the structure and emission from the transition zone of Balmer- 
dominated SNRs depend on shock velocities (ranging from v s = 300 — 10,000 km s _1 ), 
pre-shock ionization fractions (ranging from f p — 0.1 — 0.9), and temperature equilibration 
ratios (ranging from /3 = 0.1 — 1), as well as Case A and Case B conditions for the narrow 
line. 



5.1. Spatial Structure of the Shock Transition Zone 



The basic features of the shock transition zone density structure were shown in Figure [2J 
The details of broad neutral production are highly sensitive to v s , f p , and (3. We can see 
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the effect of increasing shock velocity by comparing Figure [6] to Figure [21 At high v s , the 
ionization rate is greater than the charge transfer rate, so that relatively few broad neutrals 
are produced. The densities of the broad neutral populations decrease rapidly with each 
subsequent charge transfer reaction. The helium ionization rates decrease, causing the He + 
and alpha particle production to peak further downstream. 

In Paper 2, we explored variations in ion velocity throughout the shock transition zone 
and found that, for v s > 300 km/s, Uj ~ v s /4 with negligible deviation. We confirm this 
conclusion with our multi-component models. Figure [7] shows the dimensionless ion velocity 
and temperature as a function of distance from the shock front. Results are shown for 
v s = 1000 km s _1 and v s = 4000 km s -1 , with f p = 0.5, (3 = 1. In the left panel, we 
display the percent deviation of the ion velocity from v s /A for the two shock velocities. At 
v s < 1600 km s _1 , broad neutrals are produced with smaller momentum density than the 
original ion population, leading to an increase in ion velocity by conservation of momentum 
(solid curve). For v s > 1600 km s" 1 , the opposite occurs, and the ion velocity decreases as 
the broads are produced (dotted curve). The deviation of ion velocity from v s /A is less than 
1% and has a negligible effect on calculations of the reaction rate coefficients . In the right 
panel, we plot the ion temperature profile. In this case, the broad neutrals are produced 
with slightly lower temperature than the ions, leading to slight heating of the ions. At large 
distances downstream, the presence of alpha particles increases the mean atomic mass in the 
gas to 1.27 m p (assuming a 10% helium abundance), leading to a final ion temperature that 
is slightly higher than e.; ~ 3/16(1 + (3). The maximum deviation of from its expected 
value is of order 10%, which in practice has a negligible effect on the values of the reaction 
rate coefficients. 



5.2. Ha Emissivity 

In Figure [31 we show a typical emissivity profile at relatively low shock velocity and 
ionization fraction. In the left panel of Figure El we show the effect of increasing the initial 
ionization fraction to f p = 0.9. We see that the neutral density decreases rapidly and the 
narrow line emission peaks at the shock front (left panel). Many charge transfer reactions 
occur close to the shock, pushing the centroid of the broad emission farther upstream com- 
pared to the f p = 0.1 case. Comparing the right panel of Figure M to Figure El we can see 
the effect of increasing the shock velocity. For the high velocity case, both ionization and 
charge transfer rates are decreased, shifting the centroids of both line components farther 
downstream in the transition zone. 



As discussed by 



Raymond et al 



(120071 ) and Paper 2, the spatial shift between the broad 
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and narrow line centroids potentially provides a constraint on the pre-shock ionization frac- 
tion, f p , and external density, n . In Figure [9], we plot z s h = z cb — z cn as a function of shock 
velocity v s , with n = 1 cm -3 . Note that the spatial shift scales as z oc l/n . In the left 
panel, we show how z sil depends on f p , holding (3 — 1. As the velocity increases, the charge 
transfer rate decreases relative to the ionization rate, delaying the production of broad neu- 
trals. Consequently, the centroid of £& shifts downstream and the value of z s h increases. For 
small f p , few protons initially exist to engage in charge exchange. Consequently, z c b shifts 
downstream from the shock front. In the right panel, we show how the results depend on 
the temperature equilibration parameter, f3, for fixed f p = 0.5. At shock velocities v s > 1000 
km s _1 , increasing /3 reduces the proton ionization rate for cold hydrogen. This effect causes 
the peak of broad production to shift downstream from shock front. At high velocities, 
increasing (3 tends to increase the charge transfer rate for broad neutrals over that for cold 
hydrogen. This effect causes the peak of broad production to shift closer to the shock front, 
and decreases the value of 2 s h- For more discussion on the observational significance of the 
spatial emissivity profile and the shift, see §6.21 

Another important observational diagnostic for Balmer-dominated SNRs is the ratio of 
integrated broad-to-narrow line strengths, defined as: 

f^°dz^ b (z) 

h/in= r oo , rrr- ( 36 ) 

Jo dz£ n (z) 

This ratio has a strong dependence on both v s and (3. Figure [TU] shows I\,/I n versus shock 
velocity at fixed f p = 0.5, for several values of temperature equilibration (3 — 0.1, 0.5, 1, using 
Case A (left panel) and Case B (right panel) conditions for the narrow line. For a given (3, 
the variation of the intensity ratio with velocity is the result of competition between charge 
transfer and ionization, which contribute equally at v s ~ 2000 km s _1 . At very low v s , the 
ionization rate begins to decrease precipitously while the charge transfer rate stays roughly 
constant, leading to the spike in I&/I n . At a fixed shock velocity, increasing (3 causes a 
decrease in the proton ionization rate for broad neutrals compared to cold hydrogen, leading 
to an increase in the intensity ratio. In Case B conditions (right), additional narrow emission 
due to absorption of trapped Ly/5 photons decreases the values of h/I n relative to Case A. 
As noted in SJH for most Balmer-dominated SNRs, partial line scattering of Lj(3 photons 
yields emission intermediate between Cases A and B. 

According to the picture in Paper 1, the dependence of h/In on f p is expected to be 
weak, since the number of broad neutrals produced in each population per cold hydrogen 
atom is fixed. This intuition is confirmed by our multi-component calculation. When f p 
is increased, we expect more charge transfers in the shock transition zone. However, this 
effect is balanced by higher ionization rates for both broad and cold neutral populations. 
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The net result is a negligible change in the h/I n ratio. The presence of neutral helium 
introduces a weak dependence of I b /I n on f p . In Figure [TTJ we show I b / I n as a function 
of v s for fixed (3 — 1 in Case A (left panel) and Case B (right panel) conditions. As f p 
decreases, fewer protons engage in charge transfer and ionization reactions close to the shock 
front, shifting the peak of broad neutral production downstream. In this case, the broad 
neutrals persist far enough downstream to interact with the charged helium species produced 
there, altering the ratio h/In- This effect is of order < 1% for fu e = 0.1 and remains small 
for larger helium fractions. However, it should be noted that, when the effects of partial 
Ly/3 scattering are included, variations in optical depth with f p ca n introduce a dependence 



of lb/ I n on pre-shock ionization fraction fjGhavamian et al.ll200ll ). While we have roughly 



treated Ly/3 scattering using the prescription of §H a full calculation of line scattering in the 
shock transition zone is needed to properly take into account these effects and incorporate 
the dependence of the integrated line ratio on the pre-shock ionization fraction. This is a 
source of systematic error in our calculation. 

Neglecting the f p dependence, the broad-to-narrow intensity ratio can be written as a 
function 

I b /I n = L bn {v s ,f3), (37) 

where h/In and L bn represent the measured and theoretical values of the intensity ratio, 
respectively. 



5.3. Broad Line Profiles and Observational Interpretations 

Eqs. (1321) and (|33|) allow us to model the FWHM of the broad line as a function of v s 
and /3, as shown in Fig. [5j In the extreme scenarios of Case A and Case B conditions, the 
FWHM has a weak (< 1%) dependence on the pre-shock ionization after taking a spatial 
average over the shock transition zone. Figure [12] shows examples of broad neutral velocity 
distributions [which are inputs to eq. (132]) ] for the case of face-on orientation of the shock 
as a function of line-of-sight velocity v r , in a frame of reference where the ion velocity is 
zero. Results are shown for fixed f p = 0.5 at several values of (3 — 0.01,0.1,0.5, 1 at shock 
velocities v s = 1000 km s -1 (left panel) and v s = 7000 km s _1 (right panel). At low shock 
velocities, charge transfer is extremely efficient, yielding broad neutral distribution functions 
nearly identical to that of thermal protons, with the broad neutral moving slightly slower 
than the ions. At high shock velocities, the broad neutral distributions are skewed and offset 
from the proton distribution, leading to the asymmetric profiles depicted in the right panel, 
with the broad neutrals moving considerably faster than the ions. 
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We write the FWHM relation as a function of v s and ft: 

W FWRM = W(v s ,{3), (38) 

where Wfwhm and W represent the measured and theoretical values of the line FWHM, 
respectively. For SNRs with measured Wfwhm and h/In, combining eqs. (1371) and (1551) 
self-consistently constrains v s and ft. This is accomplished by inverting W(v 8 ,ft) to yield 
v s = W /_1 (W / fwhm, ft), and calculating the root of the expression 

Un [W- l (W FWnM ,ft),ft] - h/ In = 0. (39) 

This procedure produces a pair of values (v s ,ft) for which the theoretical calculations L bn 
and W equal the measured I&/I n and W-fwhm- 

We emphasize that a self-consistent calculation is critical for accurately determining the 
values of v s and ft. In the previous literature, two bracketing values for the temperature equi- 
libration are sometimes chosen (e.g., ft = 0.1, 1), and t he FWHM relation is e mployed to give 



a range of possible values for the shock velocity (e.g.. lHeng fc McCrayl 120071 ). However, this 
procedure has a major flaw. In practice, the self-consistent quantity L bn [H /_1 (H / fwhm, ft), ft] 
will have a minimum value over the range ft £ (m e /m p , 1). If h/In for a particular observed 
shock is less than this minimum value, then no pair (v s , ft) will yield the observed values 
lb/ I n and W-fwhm- In such a case, the model breaks down and quoting a range of possible 
shock velocities for two bracketing values of ft is inappropriate. Additional physics must be 
invoked to account for the observations. Moreover, when no measurement of I b /In exists, 
the bracketing procedure may or may not yield an accurate estimate of the range of shock 
velocities. 

We use the methodology described above to self-consistently extract shock parameters 
from observations of Balmer-dominated SNRs. We summarize our results in Table [U, which 
lists, from left to right, the object, reference, Ha Hfwhm and h/I n , calculated values for v s 
and ft from Ha, Lyft H^fwhm, and calculated value for v s from Lyft. If our models do not 
yield a fit to the observations, we do not list values for the shock velocity and temperature 
equilibration ratio. Our calculations show a characteristic range of ft between 0.01 and 0.1 
for v s > 1000 km s _1 . Therefore, in the case of Lyft observations where no measurement 
of h/In exists, we report derived shock velocities for the range ft = 0.01 — 0.1; further 
observations are required to confirm the accuracy of these estimates. For SNR 0519-69.0, 
observations exist in Ha and Lyft. To calculate shock velocities from the Lyft observations, 
we use the derived value of ft from the Ha diagnostics. 



Our models successfully fit 14 measurements from seven Balmer- dominated SNRs, within 
the observational uncertainties. In the cases of the Cygnus Loop (jGhavamian et al.ll200ll ) 
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and one measurement from 0519-69.0 by lTuohy et al.1 (Il982l ). the observed ratio h/In is too 
low to be accounted for by our calculations. In addition, we are unable to fit the majority of 
measurements from the object DEM L71/0505-67.9, which we have omitted from the table. 
Below, we discuss a possible explanation for these discrepancies in our model. 

We note that our inferred shock velocities and temperature equilibration ratios are sys- 
tematically lower than in previous studies, and show an increased sensitivity to compared 
to Papers 1 and 2. The ability to sensitively probe both shock velocity and temperature 
equilibration is of interest for studies of collisionless electron heating in shocks, as described 
below. Typically, inferred values for v s are ~ 10 — 30% smaller than those quoted in Paper 
1, primarily because of the contribution of broad neutral velocities to relative speeds in fast 
neutral-ion interactions. The shock speeds in Tycho's SNR and SN1006 are smaller by about 
15% and 27%, respectively, than those reported previously. One implication is that the dis- 
tances to these SNRs derived from the shock speeds and proper mot ions are correspondin gly 



smaller. The inferred distance to SN1006 is reduced from 2.18 kpc (IWinkler et al. 



20031 ). to 



1.6 kpc. The corresponding brightness of the SN is squarely in the middle of the Type I SN 
distribution at 2.18 kpc, but is 0.7 magnitudes fainter at 1.6 kpc. The smaller distance is 
more comfortable in comparison with that of the S-M star which lies between 1.05 kp c and 
2.1 kpc and whose spectrum shows absorption by SN1006 ejecta (IBurleigh et al.1 120001) . O n 
the other hand, the ejecta are observed to expand at 7026 km s -1 f lHamilton et al.1 120071 ). 
and the requirement that this material lies within the remnant places a conservative lower 
limit to the distance of 1.6 kpc, which is just consistent with that derived here. 



Discussion 



6.1. Collisionless Shock Heating 



One of the motivations for studying Balmer-dominated SNRs is to probe the physics 
of collisionless shocks. While various mechanisms for transferring energy from proton to 
electron populations have been proposed in the literature, there is no consensus on how 
to predict the value of for a gi v en set of shock parameters in astrophysical contexts. In 
a recent paper, iGhavamian et al.1 (l2007bl ) attempted to address this problem by deriving a 
relation for temperature equilibration ratio versus shock velocity from observations of shocks 



in SNRs. They reported that their results could be fit by a curve 0(v s 



oc v a 



with 



for v s < 400 km s _1 . Given that the proton temperature scales roughly as T p oc v 2 s) this 
conclusion implies that the e lectrons are heated to a c onstant temperature, independently 
of the shock velocity. While IGhavamian et al.1 (l2007bl ) discussed a possible mechanism for 
this dependence, it has not yet been established that theoretical models can produce such 
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an effect. 

Using our new model for emission from Balmer-dominated shocks, we have calculated 
an updated (3{v s ) relation, which is displayed in Figure [T3J The solid curve depicts the 
proposed v~ 2 dependence, with (3 = 1 for v s < 400 km s _1 . For all the points shown in 
the plot, we have self-consistently fit measurements of both Wfwhm and h/I n . We fix the 
pre-shock ionization at f p = 0.5, and exclude the majority of measurements from DEM 
L71/0505-67.9, as well as one measurement each from Cygnus and 0519-69.0 which we are 
unable to account for with our models. We find from our calculations that the data are not 
well fit by a power-law relation (3{y s ) cx v~ a - if we set a = 2, the fit to the (3 curve yields 
a reduced chi-square of Xr — 62.8/12 = 5.2o For shock velocities v s < 1500 km s _1 , our 
results are fairly close to those of previous models. At higher velocities v s > 2000 km s _1 , the 
deviations are more significant. The minimum value of the temperature equilibration ratio is 
(3 ~ 0.03 at velocities v s > 1500 km s -1 . This value is greater by several orders of magnitude 
than the theoretical minimum f3 = m e /m p , but is smaller than the value predicted by some 



collisionless heating models (e.g.. ICargill fc Papadopouloslll988l ). 



6.2. Limitations to our Model and Future Work 



As demonstrated by SNR 0505-67.9, the Cygnus Loop, and 0519-69.0, an additional 
physical mechanism is needed to account for the observed Ha emission seen in some SNRs. 
One possibility is that a significant amount of the dissipated energy in the shock is transferred 
to cosmic rays, producing a pre cursor which can heat and accelerate t he upstream gas, alter - 
ing the shock jump conditions (ISmith et al.lll99ll ; iHester et al.lll994l ; ISollerman et al.ll2003l ). 
In addition, the precursor can "push" on the upst ream protons, leading to a velocity differ- 



ential between neutrals and charged species (e.g.. iBerezhko &: Ellison! Il999l ). Furthermore, 
all previous models have assumed that the kinetic distribution functions for protons and 
electrons in the shock transition zone are Maxwellian. However, recent calculations have 
shown that the proton distribution can s ignificantly depart from Maxwellian behavior, form- 
ing a distinct "pickup" ion population ([Raymond et al.ll2008l ). This effect will change the 



structure of the transition zone, broad neutral distributions, and kinetic reaction rates. 

The inclusion of a cosmic ray precursor will affect the predicted broad-to-narrow line 
strength in several ways. Broadening of the cold neutral distribution function effectively 
reduces the optical depth to Ly/3 scattering, which will decrease the conversion efficiency of 



2 Note that such a large value of \r should be interpreted in the context of 10% uncertainties associated 
with the atomic cross sections used. 
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narrow Ly/3 to Ha and increase L&„. In contrast, additional excitation of cold neutral H atoms 
in the precursor will increase narrow line emission and decrease Lf, n . While it is not obvious 
which effect is dominant, evidence for such precursor e ffects exists in the ano malously large 
widths of Ha lines in most Balmer-dominated SNRs (jSollerman et al.ll2003l ). In addition, 
H a emission f rom a spatially resolved precursor in Tycho's SNR has recently been reported 
bv lLee et all (hoOTT ). 



The spatial structure of the shock transition zone provides a way to infer the external 
density Uq and pre-shock ionization fraction f p . As noted in §5.21 the spatial emissivity pro- 
file and centroi d shift between t he na rrow and broad components have a strong dependence 
on f p and n . iRaymond et al.l (120071 ) were able to spatially resolve the total Ha emission 
from SN1006 using the ACS camera on the Hubble Spac e Telescope. iRavmond et al.l (120071 ) 
adopted the values of v s and /3 from those inferred by iGhavamian et al.l (120021 ) using ob- 
servations of the FWHM and h/I n , and then attempted to fit for the spatial structure by 
varying f p and no. Future instruments similar to ACS with high spatial resolution (e.g., 
WFC3 onboard Hubble with the narrow Ha component isolated by a narrow band filter) will 
allow for more accurate measurements of the transition zone structure and provide additional 
constraints on shock parameters. 

The models presented in our series of papers generically describe the physics of non- 
radiative shocks interacting with cold, pre-shock gas. In addition to using them to un- 
derstand nearby SNRs, we can also calculate hydrogen line emiss ion from SNRs in young , 
distant galaxies (with redshifts 3 — 5), a subject first explored by iHeng fe Sunyaevl (120081 ). 
We improve on their estimates of the luminosity ratios of Lya a nd Ly/3 to Ha, denote d 
F fy /ffn and Tr, u fi /Hry, r espectively (see Figure [Til c.f. Figure 1 of IHeng fc Sunyaevl |2 0081 ) . 



Heng fc Sunyaevl (120081 ) underestimated the broad Lya emission at high shock velocities due 
to their neglect of the mult i- component shock transition zone. Nevertheless, their conclu- 
sions remain intact: the sensitivity of T Lya / Ha and T Lyf3 / Ha to (3 by a factor ~ 2 over the 
velocity range 1000 < v s < 4000 km s _1 is a direct and unique way to measure the tem- 
perature ratio. A valuable extension to their work will incorporate models of Lya and Ly/3 
scattering, taking into account geometric effects. 

The biggest systematic uncertainties in our current model are: i) a proper treatment of 
Lyman line scattering in the shock transition zone; ii) the effect of a cosmic ray precursor 
on Lb n ; and iii) inclusion of a non-thermal population of protons in calculating the reaction 
kinetics. The incorporation of these physical effects is important for future models. Our new 
results yield substantial differences in derived shock speed (and hence inferred distances) 
for measurements of Balmer-dominated SNRs compared to previous models. While our 
results are not consistent with a power-law relation between shock velocity and electron 
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temperature, future work is needed to resolve the remaining approximations in our model 
and make a more definitive statement about shock heating of electrons. 

We thank Carles Badenes, Anatoly Spitkovsky and Eliot Quataert for useful conversa- 
tions, and Parviz Ghavamian for helping clarify several points in the observational literature. 
K.H. thanks the Institute for Advanced Study for their generous support. 



A. Shock Kinetics 



We summarize important details of our calculations of rate coefficients, broad neutral 
velocity distributions and broad line profiles. For a more detailed description of our methods, 
see Paper 1. In this work, we treat the nl atomic sub-levels separately instead of considering 
a single, summed n level as in Paper 1. We treat charge transfers, excitation and ioniza- 
tion events between elect r ons, p r otons, alpha partic le s and hydrogen atoms using the cross 



in 



(Il998h. and 



Harel et al 



Heng fc Sunyaev 



sectio ns of iBarnett et alj (119901) , iBelkic et al.l (119921 ) , Uanev fc Smith! (119931 ) , iBalanca et al. 



[ 19981 ). Fitting functions for some of these cross sections are provided 



(120081 ) 



We also consider the ionization of helium atom s and s i ngly- ionized helium by electrons 



and p r otons, using the cross sections oflPeart et al.l (119691). lAngel et al.l (119781 ) , iRudd et al. 



( 1983UShah & Gilbodv! (ll985MRinn et al.l (Il986h and lShah et al.l (Il988f ). Followin dBarnett et al 
(I1990l ). we fit the cross sections to the function 



T(x- A) 



A 



exp 



x 



(Al) 



i=i 



where the components of A = (A , A 2 , 
d are the Chebyshev polynomials: 

C (x) 
C x (x) 
d (x) 



As) are the fitting parameters, and the quantities 



(A2) 
(A3) 
(A4) 



x. 
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We define the fitting variable 



X 



In [E 2 / (E mSiX E min )} 



(A5) 



In (£ max /£ 

min ) 

where E is the relative energy between the interacting particles and E m[n and E max are the 
minimum and maximum energies for which data are available. We assume a fiducial error of 
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10% for the data. The cross sections and corresponding fits are displayed in Figure [15] and 
the fitting parameters are presented in Table [2j 

We neglect charge transfer reactions of neutral and singly-ionized helium with protons 
due to scarcity of cross sections for these processes in our velocity range of interest. While 
such interactions may affect helium line emission, we do not expect them to strongly couple 
to the Balmer radiation. In addition, for ease of computation, we approximate excitation 
and ionization of hydrogen by He + using relevant proton rate coefficients^] However, cross 
sections for this pro cess do exist in the literature and should be used in future calculations 



(IBarnett et al.lll990l ). Nevertheless, we expect corrections to these cross sections to have a 



small quantitative effect on our results. 

We make several approximations to speed up our computations. For temperature ratios 
0.01 < /3 < 0.1, the velocity width of the electron distribution is generally broader than that 
of the broad neutrals (characteristic electron velocities are greater than proton velocities by 
a factor of 4 — 14), implying that rate coefficients for interactions involving electrons are 
insensitive to changes in the broad neutral velocity distribution. We therefore approximate 
electron rate coefficients to be the same for reactions involving cold atoms and broad neutrals. 
Calculations of the FWHM for line profiles should include excitations by electrons, protons, 
singly-ionized helium and alpha particles. However, alpha particles contribute < 1% due to 
the relatively smaller density n a , and can typically be omitted. 

The reaction rate coefficients used in eqs. (J1T)- (JTT|) are defined as: 

Rx, s = J d\ J d\'f a (v)f s (V)Ava XtS (Av), (A6) 

where X denotes the interaction (ionization or charge transfer), a the atomic species (broad 
or cold neutrals), and s the interacting charged particle (electrons, protons, singly-ionized 
helium, or alpha particles). In contrast to eqs. (1531) - (1351) . the interaction cross sections a x<s 
used here represent the sum of the reactions to all nl levels for which atomic data is available 
[c.f., eq. (1271)]. Using the definition of eq. flA6j) . the quantity n a n s Rx, s gives the number of 
interactions X, between the species a and s, per unit volume. 



3 Our assumption is based on the Weizsacker- Williams approximation in w hich scattering dynamics are 

This approx- 



dominated by the charge of the impacting particle as opposed to its mass (e.g.. iJackson 
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imation is valid when the relative velocity of the collision is greater than that of an electron orbiting the 
hydrogen atom. At low velocities v 8 < 300 km s _1 , this assumption breaks down. 
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B. Typographical Errors in Papers 1 & 2 



W e point out several minor typograph ical errors in lHeng &: McCrayl (120071 ) and lHeng et al 
(b()07t ). In §5.1 of lHeng fc McCravl (|2007h. n bn (Ka) = I b (Ha)/I b (Ha) should be K bn (Yta) = 



Heng et al.l (120071 ). Rnab sh ould be iiW h n, -RHn -h. should be Rn a ,b, 



I b (Ha)/I n (Ha). In §5 of 
and Ruan should be -RHa.n- In §6.2, paragraph 2 of iHeng et al.l (l2007f l. the sentence that 
begins "Before CKR80 and HM07, ..." should be changed to "Before (CKR80 and HM07), 
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Table 1: Inferred Shock Velocities and Temperature Equilibrations from SNRs 

Object Reference Ha FWHM Ha I b /I n Ha v s Ly{3 FWHM Ly/3 







(km s 1 ) 




(km s 1 ) 




(km s" 1 ) 


(km s" 1 ) 


Cygnus 


G01 


262 ± 32 


0.59 ±0.3 


— 


— 


— 


— 


RCW 86 


G07b 


325 ± 10 


1.06 ±0.1 


367l£ 


fl £77+0.287 
v -° 1 1 -0.136 


— 


— 




G01 


562 ± 18 


1.18 ±0.03 


_ _ _ 1 Q 

690l? 31 


0.82618;8?8 





— 




G07b 


640 ± 35 


1.0 ±0.2 


673l» 


0.376_ 328 


— 


— 


0505—67.9 


S91 


580 ± 70 


>0.7 


— 


— 


— 


— 




G03 


7851^ 


o.93±g:i5 


8081™ 


0.33118;^ 


— 


— 




G03 


1055li|g 


o.88±S:iS 


1024 ± 107 


0.19418;^ 


— 


— 




G07a 


— 


— 


— 


— 


1135 ±30 


1036 ± 46* 




G07a 


— 


— 


— 


— 


1365 ± 75 


1249 ± 89* 


0548—70.4 


S91 


760 ± 140 


1.1 ±0.2 




u.^ou_ 194 


— 


— 


0519—69.0 


S91 


1300 ± 200 


0.8 ±0.2 


11781$ 


031 + - 176 

U.UOi„o ooi 


— 


— 




T82 


2800 ± 300 


0.4— 0.8 1 " 





— 


— 






G07a 










3130 ± 155 


29841™ 


Kepler 


F89 


1750 ±200 


1.1 ±0.25 


15891™ 


o n3 c i +0 - 010 

U.UOd_ 009 






Tycho 


G01 


1765 ± 110 


0.67 ±0.1 


1606+^ 


u - u ^°-0.006 








KWC87 


1800 ± 100 


1.08 ±0.16 


1634l8i 


U.UOU_ 005 








S91 


1900 ±300 


0.77 ±0.09 


i yoq+285 
1 1 OO_ 280 


47+0 .018 
u - u +' -0.011 






SN 1006 


G02 


2290 ± 80 


o.84i8;8? 


20981^ 


U.UOO_ 006 








S91 


2310 ±210 


0.73 ±0.06 


2126+ 257 


069+ - 030 








KWC87 


2600 ± 100 


0.77 ±0.08 


24261** 


n73 +0 - 125 

U.U / O_ 049 







Note: Case A conditions are assumed for the broad line. The pre-shock ionization fraction is taken to be f p = 0.5. 

*: Values are reported for the range 0.01 < fi < 0.1. 

f: We regard this as 0.6 ± 0.2 when solving for /3 and v a . 
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Table 2: Fits to Helium Species Ionization Cross Sections 





He + e" 


He+ + e" 


He + p 


He+ + p 


A, 


-78.4712 


-82.6155 


-77.0261 


-80.7740 


A 1 


-0.832236 


-0.535745 


0.596233 


1.50686 


A 2 


-1.00452 


-1.15893 


-1.37165 


-2.04982 


A, 


0.482606 


0.644513 


0.205854 


0.384422 


A A 


-0.244927 


-0.419765 


0.123038 


0.353731 


A 5 


0.121965 


0.299795 


-0.0671246 


-0.268015 


A, 


-0.0795100 


-0.204216 


-0.0135306 


-0.0107265 


A 7 


0.0537985 


0.133238 


0.0163305 


0.0101069 


A 8 


-0.0521662 


-0.0716388 


-0.00556328 


-0.110770 


E ■ 

- L 'min 


26.6 eV 


54.5 eV 


5.0 keV 


2.98 keV 




10 4 eV 


10 4 eV 


2.38 x 10 3 keV 


1.03 x 10 3 keV 
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Fig. 1. — Schematic depiction of density variation in the shock transition zone. The vertical 
axis represents dimensionless density for the different particles species, while the horizontal 
axis represents distance. The shock front is marked by the left dotted vertical line. Upstream 
of the shock, we assume that negligible interactions take place between the particles (though 
see the discussion in §2D- After protons and electrons are isotropized and decelerated at the 
shock front, ionization and charge transfer reactions deplete the cold neutrals and produce 
hot broad atoms. The transition zone terminates where all the neutral species are ionized, 
leaving a mix of protons, electrons and alpha particles. 



-30- 




Fig. 2. — Dimensionless density as a function of physical distance behind the shock front, 
for parameters v s = 1000 km s -1 , f p = 0.5, (3 = 1, n — 1 cm -3 , and fu e = 0.1. The 
left panel shows the shock transition zone structure for the neutral species; the solid and 
dash-dotted curves show the density of neutral hydrogen and helium, respectively. Charge 
transfer reactions produce broad neutrals which have undergone one (dotted curve), two 
(short-dashed) and three or more (long-dashed) charge transfer reactions. The right panel 
shows the transition zone structure for the charged species. Note that the densities of 
the protons and electrons have been scaled by a factor of 1/18. Singly- ionized helium (solid 
curve) is produced downstream from the shock and ionized to produce alpha particles (dotted 
curve). The proton (short-dashed curve) and electron (long-dashed curve) densities saturate 
when the neutral species are fully ionized. 
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12 3 

z (10 15 cm) 

Fig. 3. — Emissivity profiles for narrow and broad emission as a function of distance behind 
the shock front. Results are shown for v s = 1000 km s" 1 , f p = 0.1, (3=1, and /he = 0.1. 
The centroids of the narrow and broad profiles are indicated by the light vertical lines. 
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Fig. 4. — Symmetric, edge-wise broad neutral velocity distributions as a function of line-of- 
sight velocity v r , in a reference frame where the proton velocity is zero. Results are shown for 
fixed f p = 0.5 and several values of (5 — 0.01, 0.1, 0.5, 1. As (3 is increased, the temperatures 
of the proton and broad neutral distribution functions decrease, leading to smaller FWHM. 




Fig. 5. — Broad, edge-wise Ha line FWHM as a function of shock velocity for fixed f p = 0.5 
at several values of (3 — 0.01,0.1,0.5, 1 for Ha transitions. As (3 is increased, the predicted 
FWHM decreases (see Figure @J. 
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Fig. 6. — Same as Figure [2] except that v s = 4000 km s 1 . The density of the cold neutral 
hydrogen is scaled by 1/4. 




Fig. 7. — Percent deviation of the ion velocity from v s /4 (left panel) and the ion temperature 
(right panel) as a function of position behind the shock front for two values of the shock 
velocity: v s = 1000 km s _1 (solid curve) and v s = 4000 km s _1 (dotted curve). Results are 
shown for f p = 0.5, (3 = 1, and /n e = 0.1. 
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Fig. 8. — Same as Figure[3j except that f p = 0.9, v s = 1000 km s 1 (left panel) and f p = 0.1, 
v s = 4000 km s _1 (right panel). 
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Fig. 9. — Spatial shift in units of 10 15 cm as a function of shock velocity v s . The pre-shock 
density is set to no = 1 cm -3 ; note that z s h scales as 1/tiq. The left panel shows the variation 
of the shift with f p at fixed (5 — 1. From top to bottom, the curves show results for pre-shock 
ionization fractions f p = 0.1,0.5,0.9. The right panel shows the variation of the shift with 
(3, for fixed f p = 0.5. From bottom to top at v s ~ 2000 km s _1 , the curves show results for 
temperature equilibration ratios (3 = 0.1, 0.5, 1. 
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v„ (10 B cm/s) v s (10 8 cm/s) 

Fig. 10. — Broad to narrow intensity ratio If,/I n as a function of shock velocity v s with 
f p = 0.5, /hc = 0.1 for Case A (left panel) and Case B (right panel) conditions. From 
bottom to top, the curves show results for temperature equilibration ratios (3 = 0.1, 0.5, 1. 




v a (10 B om/s) v s (10 B cm/s) 



Fig. 11. — Broad to narrow intensity ratio as a function of shock velocity v s , with 

(3=1, fue = 0.1, Case A (left panel) and Case B (right panel) conditions. From bottom to 
top, the curves show results for pre-shock ionization fraction f p = 0.1, 0.5, 0.9. 
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Fig. 12. — Face-on broad neutral velocity distributions as a function of line-of-sight velocity 
v z , in a reference frame where the proton velocity is zero. Results are shown for fixed 
f p = 0.5, several values of f3 — 0.01,0.1,0.5,1, at shock velocities v s = 1000 km s -1 (left 
panel) and v s = 7000 km s _1 . At low shock velocity, charge transfer is extremely efficient, 
and the broad neutral distributions are very close to the thermal proton distribution, with 
the broad neutrals moving slightly slower than the ions (left panel). At high shock velocity, 
the broad neutral distribution function is skewed and offset from that of the protons, leading 
to asymmetric velocity profiles with the broad neutrals moving considerably faster than the 
ions (right panel). 
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Fig. 13. — Temperature equilibration ratio (3 as a function of shock velocity v s for Balmer- 
dominated shocks fit by our calculations. The symbol shape denotes the origin of the data 
point from SNRs in our sample. The solid curve shows the dependence (3{v s ) oc v~ 2 for a 
heating mechanism which is independent of shock velocity. Fitting this curve to the inferred 
values from our new model yields a reduced chi-square of xl — 62.8/12 = 5.2. 
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Fig. 14. — Luminosity ratios of broad Lya and Ly/3 to Ha, denoted Tiya/Ha and Tiyp/u, 
respectively. 
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Fig. 15. — Cross sections for ionization of neutral and singly-ionized helium by electrons and 
protons. For references to the atomic data, see Appendix [A] 



